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Coastal  Foredune  Evolution,  Part  2: 
Modeling  Approaches  for  Meso-Scale 

Morphologic  Evolution 

by  Margaret  L.  Palmsten1,  Katherine  L.  Brodie2, 

and  Nicholas  J.  Spore2 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  is  the  second 
of  two  CHETNs  focused  on  improving  technologies  to  forecast  coastal  foredune  evolution. 
Part  1  summarized  the  short-,  meso-,  and  multi-decadal-timescale  environmental  factors  and 
forcing  processes  that  influence  the  morphodynamic  evolution  of  coastal  foredunes  and 
described  a  number  of  observational  studies.  Part  2  summarizes  the  existing  approaches  used  to 
model  foredune  evolution  and  provides  a  new  probabilistic  modeling  framework  that  integrates 
hydrodynamic  and  aeolian  processes  to  forecast  foredune  evolution. 

INTRODUCTION:  As  discussed  in  Part  1,  coastal  foredunes  along  developed  coastlines  are 
dynamic  features  that  continually  evolve  in  response  to  changing  coastal  dynamics,  aeolian 
processes,  and  management  decisions  (e.g.,  Swann  et  al.  2015).  Predicting  foredune  evolution 
over  time  is  important  to  coastal  engineers  and  managers  because  foredunes  provide  ecosystem 
services  and  can  reduce  stonn  damages  to  coastal  infrastructure,  both  of  which  increase  the 
resiliency  of  coastal  areas.  A  recent  report  documenting  dune  management  challenges  on 
developed  coasts  (Elko  et  al.  2016)  specifically  identified  the  need  to  expand  both  observations 
and  modeling  approaches  across  meso-timescales,  which  are  more  directly  related  to  the 
operational  scale  of  coastal  dune  and  beach  management  decisions  and  programs  (e.g.,  over 
months  to  decades).  The  challenge  lies  in  both  incorporating  significant  short-term  processes  that 
affect  medium-  and  long-term  evolution  as  well  as  validating  month-to-decadal  evolution  models 
with  observations.  Part  1  of  this  two-part  series  summarized  the  significant  processes,  controlling 
factors,  and  morphologic  response  on  short-,  meso-,  and  multi-decadal-timescales;  highlighted 
significant  meso-decadal-timescale  observations;  and  concluded  with  a  summary  of  significant 
findings.  Part  2  builds  upon  this  knowledge,  reviews  existing  numerical  models  for  dune 
evolution,  and  recommends  a  modeling  framework. 

STATE-OF-THE-ART  DUNE  MODELS:  As  reviewed  in  Part  1,  dune  formation  and  erosion 
are  the  result  of  a  number  of  feedbacks  among  competing  processes.  Reflecting  the  abundance  of 
naturally  competing  processes,  numerous  dune  models  have  been  developed  to  address 
individual  aspects  of  dune  morphology  evolution.  The  present  state-of-the-art  in  dune 
morphology  modeling  is  focused  on  integrating  these  models  to  produce  a  single  model  that 
encompasses  dune  accretion  and  erosion  processes.  Within  the  discipline  of  quantitative  dune 
morphodynamic  modeling,  there  are  a  number  of  approaches  applied  including  process-based 
numerical  models,  parametric  models,  and  probabilistic  models.  Part  2  first  addresses  conceptual 
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models  of  dune  morphology  evolution,  then  parametric  models,  and  then  reviews  recent 
developments  in  integrated  approaches  to  model  dune  morphodynamics  including  process-based 
numerical  and  probabilistic  models,  as  well  as  describe  several  applications  of  the  models. 

Conceptual  Models.  Conceptual  models  for  coastal  dune  morphology  provide  the 
underpinnings  for  many  numerical  models  as  well  as  a  framework  for  conducting  field 
observations.  A  number  of  conceptual  models  describe  aspects  of  dune  geomorphology  (e.g., 
Hesp  2013;  Sallenger  2000),  whereas  Houser  et  al.  (2015)  described  a  more  holistic  model  for 
dune  morphology.  Each  of  these  conceptual  approaches  is  described  below. 

One  of  the  most  widely  applied  conceptual  models  is  the  stonn  impact  model  of  Sallenger  (2000), 
which  qualitatively  describes  the  relationship  between  hydrodynamic  forcing  and  dune  response 
and  provides  the  basis  for  the  process-based  numerical  model  XBeach  (Roelvink  et  al.  2009)  and 
the  probabilistic  model  of  Plant  and  Stockdon  (2012).  The  storm  impact  model  describes  four 
regimes  to  scale  the  effect  of  storm  influence  on  barrier  island  response.  In  the  swash  regime,  wave 
runup  does  not  exceed  the  base  of  the  dune.  In  the  collision  regime,  runup  exceeds  the  base  of  the 
dune  but  not  the  dune  crest.  In  the  overwash  regime,  runup  exceeds  the  crest  of  the  dune.  In  the 
inundation  regime,  the  surge  level  continuously  exceeds  the  elevation  of  the  barrier  island. 
Sallenger  (2000)  argued  that  the  magnitude  and  direction  of  sediment  transport  are  unique  for  each 
regime.  Stockdon  et  al.  (2007)  demonstrated  that  the  storm  impact  model  does  have  the  capability 
to  hindcast  the  dune  response  to  storm  forcing. 

Hesp  (2013)  described  the  evolution  of  transgressive  dune  systems  based  on  visual  analysis  of 
coastal  dune  fields  from  around  the  world.  He  categorized  dunes  into  those  with  high  sediment 
supply  and  little  vegetation  such  as  in  arid  climates  and  those  with  lower  sand  supply  and  more 
vegetation  in  temperate  and  tropical  climates.  He  described  the  evolution  of  the  backshore  dune 
field  including  stabilization  by  vegetation,  as  well  as  development  of  the  dune  field  following 
erosion  events,  and  the  behavior  of  parabolic  dunes. 

Houser  et  al.  (2015)  described  a  conceptual  model  for  barrier  island  equilibrium  states, 
depending  on  elevation.  In  the  model,  stonn  surge  parameters  including  frequency  of  occurrence 
and  elevation  relative  to  the  island  represent  a  tipping  point  between  banier  islands  with  high 
foredunes  and  those  with  low,  overwash-dominated  features.  The  type  and  growth  rate  of 
vegetation,  along  with  sediment  supply,  control  recovery  of  dunes  and  barrier  island  elevation. 

Parametric  Models.  Parametric  models  are  a  quantitative  alternative  to  more  simple 
conceptual  models.  There  is  presently  no  combined  parametric  model  for  dune  evolution,  but 
there  are  a  number  of  parametric  models  available  to  address  different  aspects  of  dune 
morphodynamics.  Here,  five  parametric  models  are  reviewed  that  address  different  aspects  of 
dune  erosion  and  accretion.  The  models  range  from  time-independent  scaling  relationships  to 
time-dependent  models  incorporating  simplified  physical  relationships. 

To  address  dune  erosion,  Long  et  al.  (2014)  identified  time-independent  scaling  relationships  that 
linked  the  magnitude  of  dune  response  to  storm  characteristics.  Specifically,  Long  et  al.  (2014) 
compared  the  change  in  elevation  of  the  dune  (Az)  normalized  by  the  height  of  the  dune  ( Dcrest ) 
to  the  difference  between  the  runup  height  ( Rhigh )  and  the  dune  height  ( Dcrest )  normalized  by  the 
significant  swash  elevation  ( S ): 
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This  scaling  relationship  explained  47%-73%  of  observed  dune  erosion  for  the  four  storms 
presented  in  Long  et  al.  (2014).  However,  the  scaling  was  not  consistent  among  the  individual 
sites  or  storms  investigated. 

Palmsten  and  Holman  (2012)  developed  a  parametric  time-dependent  model  for  dune  erosion 
during  the  collision  regime  drawing  from  the  wave-impact  models  of  Overton  and  Fisher  (Fisher 
et  al.  1986;  Overton  and  Fisher  1988;  Overton  et  al.  1994a,b)  and  large-scale  laboratory 
observations  of  dune  erosion.  The  model  relates  the  volume  of  sediment  eroded  from  the  dune  to 
the  elevation  of  the  extreme  wave  runup,  the  elevation  of  the  dune  base,  and  the  period  that  the 
dune  is  exposed  to  wave  runup  as  modeled  using  a  nonnal  distribution  and  the  Stockdon  et  al. 
(2006)  parameterization  for  wave  runup.  Splinter  and  Palmsten  (2012)  found  that  this  parametric 
approach  required  no  calibration  yet  accurately  tracked  the  position  of  the  dune  base  when  used 
in  hindcast  mode  on  field  data.  Palmsten  and  Holman  (2011)  were  also  able  to  relate  erosion 
during  the  collision  regime  to  moisture  content  of  the  dune  and  sediment  shear  stress. 

Three  models  for  dune  growth  through  aeolian  transport  are  discussed  next.  Baas  and  Nield  (2007) 
represented  vegetated  coastal  dune  growth  and  migration  in  the  Discrete  Ecogeomorphic  Aeolian 
Landscapes  model  (DECAL)  as  a  cellular  automaton,  which  is  an  extension  of  the  model  for 
unvegetated  dunes  developed  by  Wemer  (1995).  In  this  approach,  simple  rules  for  the  probability 
of  sediment  transport,  angle  of  repose,  and  effectiveness  of  vegetation  at  reducing  the  probability 
of  sediment  transport  are  used  to  model  foredune  morphodynamics.  Baas  and  Nield  (2007)  do  not 
make  detailed  calculations  of  airflow  or  sediment  transport.  The  DECAL  model  qualitatively 
reproduced  vegetated  dune  landscapes  (Baas  and  Nield  2007;  Baas  and  Nield  2010). 

Bauer  and  Davidson-Arnott  (2003)  described  a  modeling  framework  that  includes  simple 
sediment  transport  equations  and  schematic  beaches  to  describe  the  relationship  between  critical 
fetch,  maximum  fetch,  beach  geometry,  and  angle  of  wind  approach.  Using  this  highly 
schematized  approach,  the  sediment  transport  rate  is  a  tradeoff  between  increased  transport  rate 
into  the  dune  due  to  longer  fetches  and  decreased  transport  rate  due  to  more  sediment  moving 
parallel  to  the  dune  as  wind  angle  relative  to  the  dune  increases.  The  model  also  incorporates 
edge  effects  on  sediment  transport  rate  into  the  dune.  However,  Bauer  and  Davidson-Arnott 
(2003)  found  that  most  beaches  can  be  modeled  as  an  infinitely  long  beach  with  negligible  edge 
effects,  instead  of  simply  dominated  by  wind  angle  relative  to  fetch  distance. 

Building  on  the  work  of  Bauer  and  Davidson-Arnott  (2003),  Delgado-Fernandez  (2011) 
improved  the  modeling  framework  for  sediment  transport  based  on  fetch,  beach  geometry,  and 
wind  angle  by  incorporating  the  supply-limiting  effects  of  moisture  content,  snow,  and  ice. 
Incorporating  these  effects  resulted  in  better  hindcasts  of  total  sediment  transport  over  the  study 
area  for  the  9-month  observation  period.  Hindcast  total  transport  was  of  the  same  order  of 
magnitude  as  observations  and  represented  a  29%  improvement  over  neglecting  the  supply 
limiting  effects. 
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Process-Based  Numerical  Models.  Process-based  models  numerically  solve  the  equations 
for  hydrodynamics  and  sediment  transport  that  describe  dune  morphodynamics.  Presently,  no 
single  process-based  model  incorporates  both  dune  erosion  and  accretion,  although  there  are  two 
projects  underway  to  develop  erosion/accretion  models.  Here,  these  two  efforts  are  reviewed, 
and  a  description  is  provided  of  a  process-based  model  for  cross-shore  sediment  transport  that  is 
incorporated  into  an  operational  risk  assessment  tool. 

WindSurf  is  a  process-based  dune  morphology  model  presently  undergoing  development  at  the 
Technological  University  of  Delft/Deltares  in  conjunction  with  the  Building  with 
Nature/ZandMotor  project.  WindSurf  will  ultimately  combine  the  extreme  Beach  model  (XBeach) 
(Roelvink  et  al.  2009),  which  couples  the  hydrodynamic  effects  of  waves  and  currents  with 
morphological  changes  to  the  beach,  and  two  models  for  the  effects  of  aeolian  transport  on  the 
subaerial  beach,  Aeolis  and  the  Coastal  Dune  Model  (CDM).  The  Aeolis  model  includes  an 
advection  diffusion  equation  taking  into  account  the  effects  of  supply-limited  sediment  transport  on 
the  beach  (de  Vries  et  al.  2014).  CDM  simulates  the  role  of  vegetation  on  controlling  foredune 
crest  height  (Duran  and  Moore  2013).  The  following  paragraphs  describe  the  composition  of  these 
component  models  as  well  as  some  applications. 

XBeach  has  been  widely  applied  to  model  storm  impacts  on  beaches.  The  model  was  specifically 
designed  to  model  dune  erosion,  overwash,  and  island  breaching.  Hydrodynamics  included  in 
XBeach  are  directionally  spread  infragravity  wave  motions  in  the  surf  and  swash  zone,  including 
wave  breaking  and  wave  rollers.  Dune  erosion  is  modeled  as  an  avalanching  process  that,  while 
effective  with  tuning,  is  physically  unrealistic.  Example  applications  include  hindcasting  and 
forecasting,  morphologic  changes  during  hurricanes  in  the  United  States  (e.g.,  Lindemer  et  al. 
2010;  McCall  et  al.  2010),  Bora  storms  on  the  Adriatic  (e.g.,  Armaroli  et  al.  2013;  Harley  and 
Ciavola  2013),  east  coast  lows  in  Queensland,  Australia  (Splinter  and  Palmsten  2012),  and 
winter  storm  erosion  (e.g.,  Baart  et  al.  2015;  Vousdoukas  et  al.  2012). 

Work  is  underway  to  couple  XBeach  with  Aeolis  for  supply-limited  sediment  transport.  Sediment 
transport  in  Aeolis  is  modeled  as  one-dimensional  (ID)  advection  model  where  sediment 
concentration  is  calculated  using  a  Bagnold-type  sediment  transport  formulation.  Sediment 
concentration  is  modeled  by  wind  speed  to  the  third  power.  Erosion  is  modeled  by  a  pickup 
function,  and  deposition  depends  on  sediment  concentration.  Changes  to  the  sediment  bed  are 
limited  by  the  sediment  supply  rate  at  the  bed.  Initial  simulations  by  Aeolis  were  recently 
published  and  compared  with  linear  and  nonlinear  models  for  sediment  transport  (de  Vries  et  al. 
2014).  Early  comparisons  between  Aeolis  simulation  and  field  observations  at  ZandMotor  were 
presented  at  American  Geophysical  Union  Ocean  Sciences  Conference  in  2016  and  demonstrated 
the  importance  of  including  supply-limited  transport  in  accurately  simulating  sediment  transport. 

The  second  aeolian  transport  model  being  coupled  to  XBeach  is  CDM.  CDM  is  a  variation  on  a 
morphodynamic  model  for  desert  dunes  that  has  been  under  development  for  a  number  of  years 
(e.g.,  Duran  et  al.  2008;  Kroy  et  al.  2002;  Weng  et  al.  1991).  CDM  expands  on  the  previous 
versions  of  the  model  by  incorporating  the  effect  of  the  presence  of  water  on  the  threshold  shear 
stress  needed  for  sediment  motion  and  the  effect  of  coast  vegetation  on  threshold  shear  stress. 
Duran  and  Moore  (2013)  used  CDM  to  identify  scale  invariant  relationships  between  beach  profile 
and  maximum  foredune  height  for  incipient  and  mature  dunes.  The  relationships  agreed  with  a 
number  of  observations.  Duran  and  Moore  (2013)  also  found  in  their  model  that  dune  morphology 
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was  relatively  insensitive  to  details  of  vegetation  growth.  Using  field  data,  Duran  and  Moore 
(2013)  demonstrated  that  relative  dune  height  correlates  with  wave  height,  where  wave  height 
controls  the  distance  between  the  shoreline  and  the  vegetation  line,  hence  the  foredune  height. 

An  alternative  approach  in  development  at  the  U.S.  Army  Corps  of  Engineers  (USACE),  Engineer 
Research  and  Development  Center  (ERDC),  couples  the  CShore-C15  (Cl 5)  morphology  model 
with  a  simple  vegetation  and  dry-sand  growth  model  to  predict  dune  growth  following  overwash 
events.  C15  is  a  version  of  the  ID  CShore  morphology  model  (Kobayashi  et  al.  2008),  where 
individual  cross-shore  profiles  that  are  solved  in  ID  are  allowed  to  exchange  sediment  in  the 
alongshore  direction  driven  by  gradients  in  longshore  transport.  The  model  accounts  for  wave  and 
current  interaction,  bedload  and  suspended  loads,  and  wave-related  sediment  transport  within  the 
surf  and  swash  zones.  C15  is  computationally  efficient  and  stable  enough  for  use  in  a  stochastic 
framework  for  long  simulation  periods1.  While  a  fully  two-dimensional  version  of  CShore  termed 
C2Shore  has  also  been  developed  and  shown  to  accurately  predict  dune  recession  during  Hurricane 
Katrina  on  Ship  Island,  Mississippi  (Johnson  and  Grzegorzewski  2011),  C15  was  developed  to 
simulate  longer,  management-relevant  periods  while  retaining  sediment  transport  calculations  on 
short  timescales  (order  of  seconds).  C15  is  now  being  coupled  with  a  heuristic  dry-sand  dune 
growth  model,  a  community-based  vegetation  growth  model,  and  a  mechanism  for  enhanced 
deposition  by  aeolian  transport  over  vegetation2.  The  new  system  is  designed  to  simulate  long-term 
coastal  morphological  evolution  but  is  untested  at  present. 

One  of  the  earliest  process-based  models  for  beach  profiles  including  dunes  is  Simulating  Storm- 
Induced  Beach  Change  (SBEACH)  (Larson  et  al.  1990;  Larson  and  Kraus  1989).  SBEACH  is 
based  on  the  concept  that  within  the  surf  zone  sediment  flux  is  dependent  upon  the  imbalance 
between  actual  energy  dissipation  and  dissipation  on  an  equilibrium  beach  profile.  Linear  wave 
theory  and  the  wave  breaking  model  of  Dally  et  al.  (1985)  are  used  to  calculate  dissipation. 
Wave  setup  is  calculated  as  well.  Sediment  flux  increases  exponentially  from  the  offshore 
initiation  of  sediment  movement  through  wave  breaking  to  the  plunge  point  at  the  shoreline  and 
then  linearly  decreases  from  the  plunge  point  through  the  swash  zone  to  the  limit  of  wave  runup. 
The  sediment  transport  relations  in  SBEACH  were  derived  from  two  prototype  scale  experiments 
in  large  wave  tanks.  Inputs  to  SBEACH  include  an  initial  beach  profile,  a  description  of  beach 
characteristics  including  grain  size,  and  a  time  series  of  offshore  wave  conditions.  SBEACH 
outputs  an  updated  cross-shore  beach  profile  from  gradients  in  cross-shore  sediment  transport 
and  can  be  used  to  simulate  erosion  of  the  primary  foredune. 

SBEACH  was  incorporated  within  Beach -fx  by  Gravens  et  al.  (2007).  Beach-yi  provides  a 
framework  for  evaluating  beach  nourishment  projects  and  management  decisions  over  multi-year 
timescales  and  employs  a  Monte  Carlo  simulation  of  beach  change  using  SBEACH.  Boundary 
conditions  for  SBEACH  are  based  on  historical  environmental  forcing  for  the  study  site.  The 
Monte  Carlo  simulations  incorporate  uncertainty  in  model  parameters  and  environmental 
forcing.  A  major  advantage  of  Beach-/x  is  its  probabilistic  approach  through  Monte  Carlo 
simulation,  which  is  useful  for  decision  making.  While  Beach-/x  attempts  to  predict  potential 
dune  erosion  over  meso-timescales,  it  presently  has  no  mechanism  that  allows  for  dune  growth 
over  similar  timescales. 


1  B.  Johnson  and  J.  McNinch,  ERDC/CHL,  personal  communication,  May  2016. 

2  B.  Johnson,  ERDC/CHL,  and  C.  Piercy,  ERDC/EL,  personal  communication,  May  2016. 
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Probabilistic  Models.  Data-driven  probabilistic  models  are  an  alternative  to  deterministic 
models,  and  significant  literature  exists  regarding  their  applicability  to  dune  erosion.  There  are 
several  advantages  to  applying  a  probabilistic  model.  In  particular,  probabilistic  models  have  low 
computational  cost  and  inherently  include  an  uncertainty  estimate  for  decision-making.  Plant  and 
Stockdon  (2012)  developed  the  first  probabilistic  model  for  predicting  dune  morphodynamics. 
Their  approach  was  to  quantify  and  extend  the  conceptual  storm  impact  model  of  Sallenger 
(2000)  using  a  Bayesian  Network  (BN),  which  would  be  trained  using  two  sources  of  data.  Lidar 
observations  of  foredune  morphology  on  Santa  Rosa  Island,  FL,  before  and  after  Hurricane  Ivan 
in  2004  were  included  in  the  BN.  Tide,  surge,  wave  height,  and  wavelength  included  in  the  BN 
were  hindcast  using  Delft3D,  and  wave  runup  at  the  shore  was  calculated  using  the 
parameterization  of  Stockdon  et  al.  (2006).  In  addition  to  dune-crest  elevation,  dune-base 
elevation,  dune-crest  elevation  change,  surge  height,  and  runup  height  included  in  the  Sallenger 
(2000)  stonn  impact  model,  the  BN  included  dune  width,  beach  width,  dune  position  change,  and 
shoreline  change.  Using  this  approach,  Plant  and  Stockdon  (2012)  accurately  hindcasted  dune 
change  for  Hurricane  Ivan.  Palmsten  et  al.  (2014)  adapted  the  model  of  Plant  and  Stockdon 
(2012)  for  use  on  the  Gold  Coast  of  Queensland,  AU,  to  hindcast  dune  retreat  distance.  Palmsten 
et  al.  (2013)  demonstrated  that  a  Bayesian  network  developed  at  one  geographic  location  can  be 
applied  at  a  second  location,  given  similar  morphology  and  similar  hydrodynamic  forcing. 

PROBABILISTIC  DUNE  EVOLUTION  MODEL:  The  authors  propose  a  probabilistic  model 
for  dune  evolution  combining  the  approach  of  Plant  and  Stockdon  (2012)  to  model  dune  erosion 
due  to  the  impacts  of  storm  surge  and  waves  with  the  framework  described  by  Delgado- 
Femandez  (2011)  to  model  dune  accretion  due  to  aeolian  processes.  The  framework  described  by 
Delgado-Femandez  (2011)  will  be  incorporated  into  the  BN  design  developed  by  Plant  and 
Stockdon  (2012)  (Figure  1).  The  new  BN  dune  model  will  incorporate  morphological  variables 
including  dune-crest  elevation,  dune-base  elevation,  dune  width,  and  beach  width. 
Hydrodynamic  variables  will  include  surge  height,  runup  height,  and  moisture  content.  Aeolian 
variables  will  include  wind  speed,  wind  direction,  fetch  distance  (which  depends  upon  the 
distance  between  the  shoreline  and  the  vegetation  line),  and  critical  fetch  distance.  Variables 
describing  the  duration  of  hydrodynamic  and  aeolian  forcing  events  will  be  included  in  the 
network,  and  the  network  will  be  forced  with  time-dependent  forcing  for  hindcasts  and  forecasts. 

In  this  new  BN  dune  model,  variables  are  represented  as  probability  distributions  of  prior 
observations.  Arrows  linking  the  variables  (Figure  1)  represent  the  joint  probabilities  between 
variables.  The  joint  probability  distributions  between  variables  will  be  learned  using  the 
Expectation-Maximization  (EM)  algorithm,  a  widely  used  machine  learning  technique  (e.g.,  Moon 
1996).  The  necessary  data  will  be  obtained  to  train  the  network  with  a  nearly  30-year  record  of 
variables  from  hydrodynamic,  morphologic,  meteorological,  and  Argus  video  imagery  data 
collected  at  the  ERDC  Coastal  and  Hydraulics  Laboratory  (CHL),  Field  Research  Facility  (FRF)  in 
Duck,  NC.  This  approach  will  be  used  to  evaluate  the  sensitivity  of  dune  erosion  to  input  variables 
and  to  test  the  hypothesis  that  dune  building  at  the  FRF  occurs  during  windy,  low-moisture 
conditions  when  wind  direction  is  oblique  to  the  dune  line.  Once  sensitivity  is  analyzed  for  dunes 
at  the  FRF,  the  data  will  be  incorporated  from  other  sites  to  test  the  applicability  of  the 
relationships  that  are  identified.  The  new  BN  approach  described  here  could  integrate  with  existing 
models,  which  can  provide  time-dependent  forcing  to  the  network.  Additionally,  the  new  BN  dune 
model  can  be  integrated  with  Beach-/x  to  be  used  as  a  tool  for  risk  analysis  and  decision  making. 
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Sallenger  (2000)  morphologic 
Sallenger  (2000)  hydrodynamic 
Plant  and  Stockdon  (2012)  morphologic 
Delgado-Fernandez  (2011)  morphologic 
Delgado-Fernandez  (2011)  hydrodynamic 
Delgado-Fernandez  (2011)  Aeolian 
Wilson  et  al.  (2015)  hydrodynamic 
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Figure  1 .  Schematic  diagram  of  the  new  Bayesian  network  for  dune  evolution.  Boxes  represent  probability 
distributions  of  variables,  and  black  arrows  represent  the  joint  probabilities  between  variables. 
Yellow  boxes  are  morphological  variables,  green  boxes  are  aeolian  variables,  and  blue  boxes 
are  hydrodynamic  variables.  Boxes  outlined  in  solid  black  are  variables  included  in  the  Sallenger 
(2000)  storm  impact  model.  Boxes  outlined  in  the  dashed  black  line  represent  the  extended 
Bayesian  Network  model  of  Plant  and  Stockdon  (2012).  Boxes  without  outlines  represent 
variables  in  the  Delgado-Fernandez  (2011)  framework.  The  box  outlined  in  the  black  dotted  line 
represents  a  hydrodynamic  variable  from  the  Bayesian  network  of  Wilson  et  al.  (2015). 


SUMMARY  OF  FINDINGS:  Developing  a  predictive  capability  for  foredune  evolution  that 
incorporates  dune  growth  and  erosion  processes  is  important  to  coastal  engineers  and  managers 
because  understanding  the  natural  ability  of  foredunes  to  rebuild  post-storm  can  help  improve  the 
management  of  dunes  and,  ultimately,  the  ability  to  quantify  the  resiliency  of  coastal  areas.  The 
probabilistic  modeling  framework  developed  in  this  technical  note  is  based  on  observations  of 
foredune  morphodynamics  (discussed  in  Part  1)  and  existing  models.  Significant  findings  from 
Part  2  are  summarized  below. 

•  Foredune  morphology  depends  on  the  feedback  between  hydrodynamic  and  aeolian 
processes,  yet  no  single  model  presently  exists  that  combines  both  types  of  processes. 

•  The  relationships  between  beach  morphologic  characteristics  including  elevation  of  the 
dune  base  and  dune  crest,  dune  width,  beach  width,  and  hydrodynamic  forcing  including 
surge  and  wave  runup  are  the  basis  for  most  conceptual,  parametric,  numerical,  and 
probabilistic  models  for  dune  erosion. 

•  The  relationships  between  beach  morphologic  characteristics  including  dune  elevation, 
beach  width,  and  grain  size;  aeolian  characteristics  including  wind  speed,  wind  direction, 
and  fetch  distance;  and  hydrodynamic  characteristics  including  moisture  content  and 
precipitation  are  the  basis  for  most  conceptual,  parametric,  numerical,  and  probabilistic 
models  for  dune  accretion. 
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•  Incorporating  the  effects  of  time-dependent  changes  to  morphology  is  an  essential  part  of 
all  models  for  foredune  evolution. 

•  A  new  probabilistic  modeling  framework  was  designed  with  the  aim  of  developing  a 
predictive  capability  for  foredune  evolution  over  short-,  meso-,  and  multi-decadal- 
timescales  due  to  both  hydrodynamic  and  aeolian  processes.  This  model  will  be  initially 
tested  and  evaluated  using  30  years  of  data  collected  at  CHL’s  FRF  in  Duck,  NC.  Finally, 
there  will  be  an  evaluation  of  the  sensitivity  of  the  model  to  each  variable. 

ADDITIONAL  INFORMATION:  For  additional  information,  contact  Katherine  L.  Brodie, 
Coastal  Observation  and  Analysis  Branch,  Coastal  and  Hydraulics  Laboratory,  1261  Duck  Rd, 
Duck,  NC,  27949  at  252-261-6840  x233  or  email:  Katherine. L.  Brodie(cv,usace.army. mil.  This 
research  was  initiated  and  funded  by  the  USAGE  Coastal  Inlets  Research  Program  in  response  to 
Statements  of  Need  2014-N-ll,  “Guidance  and  Tools  Needed  to  Calculate  Integrated  Coastal 
Resilience,”  and  2015-N-ll,  “Develop  Guidance  for  Coastal  Resilience.”  This  CHETN  should 
be  cited  as  follows: 

Palmsten,  M.  L.,  K.  L.  Brodie,  and  N.  J.  Spore.  2017.  Coastal  foredune  evolution, 

Part  2:  Modeling  approaches  for  meso-scale  morphologic  evolution  .  ERDC/CHL 
CHETN-II-57.  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development 
Center,  http://dx.doi.ors/ 10.2 1079/1 1681/21627. 
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